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ABSTRACT 


An  Algor it ha  for  position  determination  using  Loran-C  trip¬ 
lets  and  an  implementation  of  tha  algor ithm  in  the  BASIC 
language  of  the  Commodore  2001  microcomputer  are  presented. 
Modifications  of  the  geodesic  direct  solution  and  reverse 
solution  algorithms  of  p.  D.  Thomas  are  also  presented) 
these  modifications  eliminate  all  of  the  special  position¬ 
ing  and  quadrant  determination  requirements  of  the  original 
algorithms. 


A.  introduction 

Loran-C  position  determination  programs  for  the  Hewlett- 
Packard  HP-67  programmable  calculator  for  use  aboard  u.  S. 

Navy  patrol  aircraft  (P-3)  have  been  prepared  and  are  reported 
elsewhere  [Ref.  11.  During  the  development  and  testing  of 
the  HP-67  programs  the  need  for  a  high  accuracy  program  for 
cross  checking  purposes  arose.  The  program  presented  here  was 
developed  as  a  refinement  of  the  HP-67  algorithm. 

The  basic  methodology  of  this  program  is  contained  in 
the  formulas  of  Paul  D.  Thomas  [Ref.  2]  for  the  direct  and 
reverse  problems  of  oblate  spheroidal  geodesy.  The  direct 
problem  is  to  determine  the  latitude  and  longitude  of  a  point 
P2  when  the  latitude  and  longitude  of  a  point  P^  as  well  as 
the  azimuth  and  distance  of  P2  from  P^  are  known.  The 
reverse  (or  inverse)  problem  is  to  determine  the  distance, 
forward  azimuth,  and  backward  azimuth  between  two  points,  P^ 
and  P2,  when  the  latitude  and  longitude  of  both  P^  and  P2 
are  known.  In  both  problems  P^  and  P2  are  on  the  surface  of 
a  spheroid.  Unfortunately,  Thomas'  algorithms  contain  several 
restrictions  and  numerous  subcases  that  are  too  awkward  to  be 
programmed  on  the  HP-67.  For  this  reason  portions  of  these 
algorithms  were  rewritten  so  as  to  eliminate  all  of  the  original 
restrictions  and  subcases.  These  modified  algorithms  for  the 
direct  and  reverse  problems  are  presented  here  with  notations 
about  the  modifications  made. 


B.  Loran-C  Fixing  Algorithms 

The  principles  of  Loran  lines  of  position  (LOP's)  and 
fixing  are  adequately  covered  in  Reference  3  and  will  not  be 
repeated  here. 

The  basic  Loran-c  equation  [Ref.  4]  can  be  written  as 

T  =  [Ts  +  p(Ts)]  -  [TM+p(TM)]  +  [TB  +  p(TB)]  +  6  (1) 

where 

T  is  the  "indicated  time  difference"  in  microseconds/ 

TM  is  the  distance,  in  microseconds/  from  the  master 
to  the  receiver/ 

Tg  is  the  distance,  in  microseconds,  from  the  slave 
to  the  receiver, 

TB  is  the  distance,  in  microseconds,  between  the 
master  and  the  slave, 

6  is  the  assigned  coding  delay,  in  microseconds, 
and  p(T)  is  the  secondary  phase  correction,  in  micro¬ 
seconds,  for  a  surface  water  path  of  length  T. 

The  quantity 

At  =  [TB+p(TB)]  +  5 

is  a  constant  for  each  master/slave  pair.  The  following 
World  Geodetic  System  1972  (WGS  72)  values  have  been  adopted 
for  Loran-C  navigation  [Ref.  4]  : 

Vq  =  299792458  meters/second  is  the  velocity  of 
light  in  free  space, 

n  =  1.000338  is  the  index  of  refraction  of  the 

surface  of  the  earth  for  standard  atmosphere 
and  100  kHz  electromagnetic  waves. 
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ae  >  6378135.000  meters  is  the  equatorial  radius 
of  the  earth 

and  f  ■  1/298.26  is  the  flattening  factor  (l-b/ae, 

where  b  is  the  polar  radius)  of  the  earth. 

With  these  parameters  the  secondary  phase  correction  for 
an  all  seawater  path  has  been  taken  to  be  of  the  form 

P (T)  =  aQ/T  +  ax  +  a2T  (2) 

where  T  is  in  microseconds  and, 

1.  For  T>537  ysec: 

aQ  =  129.04398, 
ax  =  -0.40758, 
and  a2  3  0.00064576438. 

2.  For  T<537  nsec: 

aQ  =  2.7412979, 
a1  =  -0.011402, 
and  a2  =  0.00032774624. 

If  one  uses  a  spherical  approximation  to  the  earth's 
surface  (f  =  0) ,  then  a  spherical  hyperbola  can  be  represented 
by  the  equation  [Ref.  3,  page  175] 

tan  r  =  2  ,  c<?s  2a  -  cos  2S, - —  (3) 

sin  2c  cos  a)  +  £  sin  2a 

where  the  origin  of  the  coordinate  system  is  at  the  prime 
focus  of  the  spherical  hyperbola,  2c  is  the  spherical  arc 
joining  the  foci,  2a  is  a  constant  for  any  one  hyperbola, 
and  r  and  co  are  the  spherical  coordinates  of  a  point  on  the 
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hyperbola.  If  the  base  line  of  the  coordinate  system  is  the 
arc  joining  the  foci  then  u>  is  the  spherical  polar  angle  from 
the  baseline  to  a  point  P  on  the  spherical  hyperbola  and  r  is 
spherical  polar  distance  (or  arc)  from  the  prime  focus  to  P. 

Using  the  Loran  system  we  take  ?  =  +1  if  the  prime  focus  is 
at  a  master  station  and  we  take  C  *  -1  if  the  prime  focus  is 
at  a  slave  station. 

If  we  take  v  =  vg/n  to  be  the  velocity  of  100kHz 
electromagnetic  radiation  at  the  earth's  surface  then,  for  a 
spherical  earth,  we  can  relate  the  parameters  in  Equations  1 
and  3  as  follows: 

2c  =  vTB/ae, 

and  2a  =  v(Tg  -  TM)/ae. 

Using  the  spherical  approximation  for  now,  we  see  that  2c  is 
known  for  any  Loran  pair.  The  "indicated  time  delay"  T  is 
measured  by  the  receiver  at  point  P,  and  to  determine  a  hyper¬ 
bolic  line  of  position  we  must  determine  2a,  but  Tg  -  TM 
cannot  be  computed  from  Equations  1  and  2.  If  a^  were  zero 
in  Equation  2,  then  it  would  be  possible  to  determine  Tg  -  TM 
uniquely.  As  a  first  approximation  we  use  the  following 
parameters  in  Equation  2: 

a0  »  o, 

a±  =  -0.321, 

and  a 2  ~  0.000635. 

These  values  have  been  obtained  by  setting  aQ  *  0  and  determining 


4 


r 


\ 


£ 

i 

■X 


and  a2  by  linear  regression  of  the  T  >  537  values  over 
the  interval  of  1000  <  T  <  8000.  This  approximation  is  quite 
good  (within  0.03  ys)  for  distances  up  to  10,000  microseconds 
where  small  changes  in  the  LOP's  can  cause  large  position 
errors.  At  short  distances  the  error  increases  from  0.05  ys  at 
1000  ys  to  0.58  ys  at  10  ys;  although  these  errors  are  large  for 
small  distances,  the  LOP's  are  not  as  sensitive  to  these 
changes  as  they  would  be  at  large  distances.  These  errors  are 
well  within  the  4  n.mi  (24  ys)  accuracy  that  was  requested  by 
COMPATWINGSPAC .  When  this  approximation  is  substituted  into 
Equation  1,  we  obtain 


[Ts  +  ax  +  a2Ts]  -  [Tm  +  ax  +  a^l  =  T  -  At, 


or 


Ts  -  tm  *  (T  -  At ) / ( 1  +  a2) 


(4) 


and  hence  2a  =  v(Tg-TM)/ae  determined  for  use  in  the 
spherical  approximation. 

Consider  a  Loran-C  triplet  with  the  master  stations  colocated. 
Let  and  E2  denote  the  azimuth  angles  of  slave  1  (S^)  and 
slave  2  (S2) ,  respectively,  measured  from  North  toward  the  East 
from  the  master  stations  (M)  (see  Figure  1).  Further,  let  a 
and  r  be  the  azimuth  and  spherical  polar  arc  (distance)  of  the 
receiver  (R)  from  M.  For  this  geometry,  Equation  3  can  be 
written  in  the  form 


tan  r^ 


Bi 


Ci 


cos 


(a  - 


*  \ 


(5) 
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where 


At  =  sin  2a ^ , 

*  cos  2a  ^  -  cos  2c^, 
and  C.  *  sin  2c^ 

for  the  Loran  pair,  i  *  1,2.  Since  ri  *  r2  *  r' 
we  can  eliminate  tan  r  between  the  two  equations.  The 
resulting  equation  can  be  rewritten  as 


C 

cos  a  +  S  sin  ot 

*  K, 

where 

C 

s  6^2  COS  ^  * 

B2C1 

cos  C., 

S 

*  Blc2  sin  ^2  “ 

B2Ci 

sin 

and 

K 

*  ®2^1  ~  ®l^2* 

If  we  define  p  >  0  and  Y  by  the  equations 

p  cos  y  -  C , 
and  psin  y  -  S, 

then 


p  */c2  +  S2, 

and  Y  *  qatn  (S,C)  . 

Here  the  function  qatn(y,x)  is  the  arctangent  of  y/x  adjusted 
for  the  proper  quadrant  according  to  the  signs  of  x  and  y. 


A  compact  form  of  this  function  is 


-1  J 

qatn(y,x)  =  tan  - o -  +  irt(x  <  0?) 

x  +  10  t (x  *  0?) 

where  t(z)  =  1  when  z  is  true 

and  t  (z)  =  0  when  z  is  false. 


When  convenient  we  will  use  the  notation  qatn(y/x)  inter¬ 
changeably  with  qatn(y,x).  Now  we  can  substitute  Eq. (7)  into 
Eq.  (6)  and  solve  for 

a  =  y  ±  cos  ^ ( K/ p )  ( 8 ) 


to  obtain  the  azimuth  angle  a  of  the  two  points  of  inter¬ 
section  of  the  spherical  hyperbolic  LOP's.  Finally  we  can 
obtain  a  value  for  r  by  substituting  each  a  into  either 
Eq.  (5) .  We  find  that 


r 


qatn 


B  • 


Ci  cos  (a-£j_)  +  Aj 


for  i 


1  or  2. 


To  obtain  the  first  solution,  r  and  a  are  entered  into 
the  "direct"  solution  algorithm;  latitude  $  and  longitude  X 
are  the  outputs. 

The  solution  is  improved  using  a  two-variable  Newton- 
Raphson  search  for  the  zeroes  of  the  functions 


f.(4>,A)  =  f.  =  (dgi  -  dM)/ae  -  2at  for  i  =  1,2, 

where  dM,  dgl,  and  dg2  are  the  distances  from  R  to  M, 
S^,  and  S2,  respectively,  and  are  computed  using  the 
"inverse"  solution  algorithm.  The  flow  of  the  improvement 
algorithm  follows: 
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1.  input  the  latest  approximation  to  the  receiver  latitude  4> 


and  longitude  X. 

2.  Compute  d^,  dg^  and  dgj  using  the  "inverse"  solution 
algorithm. 

3.  Compute  TM  *  dM/v  and  Tg^  *  dg^/v  ^or  *  *  1#2. 

4.  Compute  Tg.  -  TM  =  (T-At).  -  p(Tgi)  +  p(TM)  and 
2a j  =  v (Tsi-TM) /ae  for  i  «  1,2. 

5.  Compute  f  ±  X)  *  (dsi“dM)/,ae  -  2aA  for  i  =  1,2. 

6 .  Compute 

Sfj  _  f  ±  ( <i>  +  A*,  A)  -  f  i  (  4>,  X) 

7 4  A4 

and 

£i  _  ttU.x  )  -  tjjVX)  £or  l  .  12 

3X  AX 


-4 

(AX  =  A<j>  =  10  radians  is  used  in  the  BASIC  program.) 

7.  Replace  $  by 


and  replace  X  by 


where 


9fl 

3f  2 

8fl 

3f  2 

34> 

3  X 

3  X 

34 

8.  Repeat  from  Step  1  until  4  and  x  are  stationary. 
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C .  The  Direct  Solution  Algorithm 

This  direct  solution  algorithm  is  a  modification  of  the 
second  order  in  flattening  (f)  algorithm  given  by  Thomas 
[Ref.  2,  pp.  7-8].  Thomas1  notation  has  been  followed  as 
closely  as  possible  for  ease  of  comparison  of  the  algorithms. 
The  qatn  function  is  defined  in  the  previous  section.  West 
longitudes  and  south  latitudes  are  negative.  We  are  given  the 
point  on  the  spheroid,  where  <t>lf  \ ^  are  the  geodetic 

latitude  and  longitude  (geographic  coordinates);  the  forward 
azimuth  a12  and  distance  S  to  a  second  point  ( ^2 #  '  an<* 

from  these  we  are  to  find  the  geographic  coordinates  <J>2,X2  and 
the  back  azimuth  The  9i-ven  quantities  are  <J>^,  a^2 

and  S.  No  assumptions  about  the  relative  location  of  and  P2 
are  required.  The  modified  direct  solution  algorithm  is: 

»  tan-1[(l-f)  tan  M  =  cos  6j^  sin  a12, 

9 

N  =  cos  0^  cos  <*12'  c^  =  fM,  c2  =  f  ( 1-M  )/4 , 

D  =  (1  -  c2) (1  -  c2  -  CiM) ,  P  =  c2 [1  +  (l/2)c1M]/D  , 

=  qatn(N,  sin  0.^)  ,  d  =  S/(agD), 
u  =  2(o1-d),  W  =  1-2P  cos  u,  V  =  cos  (u+d) , 

X  =  c22  sin  d  cos  d  (2V2-1) ,  Y  =  2PVW  sin  d, 

Ao  =  d  +  X  -  Y,  Zo  -  2o1  -  A 0, 

cu.  =  qatn  [-M,  -  (N  cos  Act  -  sin  6^^  sin  Ao)  ]  , 
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K  «  (1-f)  [M2  +  (N  cos  Ac  -  sin  0^^  sin  Ao)2]1^2 
4»2  *  tan-1  [(sin  0^  cos  ha  ♦  N  sin  Ao)/K]  , 

An  *  qatn(sin  Ao  sin  ai2'cos  0^  cos  ha  -  sin  6^  sin  Ao  cos  a^2)  » 
H  *  c1(1-c2)Ao  -  c^Cj  sin  Ao  cos  Eo, 

AX  ■  An  -  H,  ^2  *  ^1  +  ^ 

In  addition  to  the  introduction  of  the  qatn  function  for 
proper  quadrant  determination  the  following  changes  have  been 
made  to  the  original  algorithm: 


t 

i 

<■ 

i 

I 


§ 

\ 


4. 

i 


1.  0q  is  no  longer  computed  from  the  equation  M  *  cos  6q. 
Since  P^  is  no  longer  required  to  be  westerly  of  P2, 
<*^2  is  no  longer  required  to  lie  in  the  interval 

[0,  180°].  Consequently  M,  which  was  originally 
required  to  be  positive,  can  be  negative.  With  this 
change  0q  was  no  longer  properly  determined  from 

M  «  COS  8g. 

2.  With  the  elimination  of  6g  the  determination  of 

from  the  equation  cos  *  cos  6^/sin  6g  became 
impossible.  The  following  figure  shows  the  spherical 
triangle  involving  and  0Q. 
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Using  the  "4-parts  formula"  [Ref.  5,  pg.  121  we  can 
determine  directly  from  the  equation 

«  qatn(cos  6^  cos  ai2'  sin  9^) 

3.  The  original  equation  for  4>2, 

tan  <J>2  *  -(sin  9^  cos  Ao  +  n  sin  Ao)  sin  a2i/d- f)M 

becomes  indeterminate  when  *  0°  or  180°  since 
both  M  and  sin  become  zero  simultaneously. 

From  the  equation  for  a2i  one  can  determine  that 

sin  a21  *  -M/[M2  +  (N  cos  Aa  -  sin  9^  sin  Ao)2]V2. 

Using  this  equation  to  eliminate  sin  <*21  in  the  equation 
for  tan  <t>2  allows  the  M  in  the  numerator  and  denominator 
to  cancel  thus  eliminating  the  indeterminancy  in  the 
equation  for  <t>2. 

With  these  changes  the  restrictions  that  be  West  of  P2» 
and  that  if  Px  and  P2  both  have  negative  latitude  the  symmetric 
positive  latitude  problem  be  solved,  have  been  eliminated. 

Also  eliminated  are  rules  for  the  determination  of  the  quadrants 
of  and  An. 

D .  The  Reverse  (Inverse)  Solution  Algorithm 

This  reverse  solution  algorithm  is  a  modification  of  the 
second  order  in  flattening  (f)  algorithm  given  by  Thomas  [Ref.  2, 
pp.  8-10].  Thomas'  notation  has  been  followed  as  closely  as 
possible  for  ease  of  comparison  of  the  algorithms.  The  qatn 
function  is  defined  in  Section  B.  West  longitudes  (X)  and 
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South  latitudes  (<j>)  are  negative.  We  are  given  the  points 

P1(,,,1'X1*'  p2  <  ^2 #  X2^  on  the  ®Pheroid  and  are  to  find  the 
distance  S  between  the  points  and  the  forward  and  back  azimuths, 

a12  and  a21*  G*ven  quantities  are  <|>3,  <j>2  and  X2-  No 

assumptions  about  the  relative  location  of  and  P2  are 
required.  The  modified  reverse  solution  algorithm  is: 

■  tan”^I(l-f)  tan  <t>^ ]  r  i  ■  1,2, 

6m  =  (0i+e2)/2'  A9m  *  (62“ei>/2'  AA  «  X2  -  X1# 

AX  «  AX/2,  H  =  cos2  A0  -sin2  0  =  cos2  0  -  sin2  A0  =  cos  0.  cos  0, 

m  m  m  m  m  l  * 

L  ■  sin2  A0m  +  H  sin2  AXm  «  sin2  (d/2), 

1-L  ■  cos2(d/2),  d  «  cos-1(l-2L), 

U  »  2  sin2  0m  cos2  A0  /(1-L),  V  *  2  sin2  A0„  cos2  0„/L, 
ro  m  ro  m 

X  »  U  +  V,  Y  •  U  -  V,  T  »  d/sin  d,  D  »  4T2, 

E  *  2  cos  d,  A  *  DE,  B  =  2D,  C  =  T  -  (A-E)/2 

nx  *  X (A+CX) ,  n2  *  Y (B+EY) ,  n3  *  DXY , 

63d  =  f  (TX-Y)/4,  <$2d  -  f2(n1-n2+n3)/64, 

S  =  ae (T-61d+62d)  sin  d,  F  *  2Y  -  E(4-x) 

M  *  32T  -  ( 20T-A) X  -  (B+4)Y, 

G  «  f T/2  +  f 2M/64 ,  Q =  -(FG  tan  AX)/4, 
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A^'  -  (  AX  ♦  Q)/2, 


tj^  ■  q»tn(-sin  A6m  cos  AA^,  cos  ^  sin  A^')» 
t2  ■  qatn(cos  A^  cos  A)^',  sin  ^  sin  A^')# 

°12  "  *1  +  t2>  "  *1  "  t2 

The  only  changes  made  to  the  original  algorithm  are  the 
computation  of  t^  and  t2  and  the  determination  of  and 

as  the  sum  and  difference  of  t^  and  t2.  Minor  though  these 
changes  may  seem,  they  have  eliminated  the  requirement  that 
is  West  of  P2  and  they  have  eliminated  four  cases  for  the 
quadrant  determination  of  and  a^.  Quadrant  determination 
is  still  required,  but  it  is  implemented  with  the  qatn 
function;  the  qatn  function  is  available  on  the  better  handheld 
calculators  in  the  form  of  the  rectangular-to-polar  function. 

E.  Program  Accuracy 

Several  programs  were  written  for  the  Commodore  2001 
computer.  One  program  implemented  only  the  direct  and  reverse 
solution  algorithms  presented  in  the  previous  two  sections. 
Initially  these  algorithms  were  validated  against  the  long 
line  computations  with  the  Clarke  1866  spheroid  model  presented 
in  Reference  2.  Later  the  reverse  solution  algorithm  was 
compared  with  the  WGS  72  data  for  40  Loran-C  pairs  contained  in 
Reference  4  (data  for  the  9930  stations  were  not  included) . 

In  these  comparisons  the  algorithm  error  for  the  baseline 
distances  was  0.16  meters  high  for  the  7980W  pair  and  0.15 
meters  low  for  the  7980Y  pair.  For  the  remaining  stations 
all  of  the  baseline  errors  were  high,  but  the  largest  error 
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was  only  0.08  meters.  It  should  be  mentioned  that  the  Micro¬ 
soft  BASIC  for  the  Commodore  2001  carries  10  digits  of  floating 
point  internally,  but  will  display  only  a  maximum  of  nine 
digits;  the  baselines  reported  in  Reference  4  are  printed  to 
the  nearest  0.0001  meters.  The  computed  one  way  baseline  time 
plus  secondary  phase  correction  was  also  computed  and  compared 
to  the  40  Loran-C  pairs  data.  In  three  cases  the  algorithm 
was  high  by  0.01  us,  in  three  cases  the  algorithm  was  low  by 
0.01  us,  and  there  was  no  error  in  the  remaining  cases. 

Another  program  was  used  to  generate  the  indicated  time 
delay  between  Loran-C  pairs  and  test  positions  for  a  receiver. 
Combining  the  time  delay  for  Loran-c  pairs  into  allowable 
triplets  enabled  the  position  fixing  program  in  the  Appendix 
to  be  evaluated.  Typical  results  are  presented  in  Table  I. 

In  all  cases  the  Commodore  2001,  running  at  BASIC  interpreter 
speed  required  8  seconds  from  entry  of  the  time  delay  input  to 
the  display  of  the  first  pair  of  solutions.  The  improvement 
algorithm  will  improve  only  one  user  selected  solution.  Each 
iteration  requires  14.4  seconds,  and  each  improvement  is  dis¬ 
played  so  that  the  user  can  monitor  the  improvement  process. 
Usually  the  solution  has  stabilized  after  two  iterations,  but 
occasionally  a  third  iteration  is  generated.  In  Table  I  the  worst 
fix  took  six  iterations  for  the  triplet  comprised  of  pairs  9930W 
and  9930X  where  the  radial  position  error  is  0.91  n.mi.  In  fair¬ 
ness  to  the  algorithm  it  should  be  noted  that  the  three  stations 
are  776  n.mi,  2130  n.mi,  and  2514  n.mi  from  the  receiver  and 
are  all  within  6°  of  azimuth  as  viewed  from  the  receiver. 


Table  I 


Receiver  Location:  52°N,  35°W 


Station 

Pair 

Computed 

Delay 

Station 

Pair 

Computed 

Delay 

7970W 

32845.12 

7930W 

16673.77 

7970X 

18889.03 

7930X 

30899.81 

9930W 

16067.09 

7930Z 

49369.30 

9930X 

28012.61 

Triad 

First 

Solution 

Improved 

Solution 

7970W/7970X 

52°00 ' 
34°59 ' 

14"N , 
34"W 

52°00' 

34°59' 

00" . 14N , 

59" . 66W 

2 

9930W/9930X 

51^56 ' 
35°21 ' 

13"N, 

23"W 

51°59* 

35°01' 

45".25N, 

25" .37W 

6 

7930W/7930X 

52°00' 
35°00 ' 

01"N , 
01"W 

51°59 ' 
35°00 ' 

59" . 69N, 

00" . 02W 

2 

7930W/7930Z 

52°00' 
35°00 ' 

01"N, 

01"W 

52°00 ' 
35°00 ' 

00" . 03N, 

00" . 07W 

2 

7930X/7930Z 

52°00' 
35  00' 

01"N » 
01"W 

52°00 ' 
34°59 ' 

00" . 00N , 

59" .96W 

2 

7930X/7970W 

52°00' 
34°59 ' 

04"N, 

57"W 

52°00' 

34°59' 

00" . 01N , 

59" . 96W 

2 

7930Z/9930X 


51°59'52"N, 
34°59 ' 26"W 


51°59'59".60N 
34°59 1 58" . 27W 
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APPENDIX.  BASIC  Program  Listing 

This  appendix  contains  the  listing  of  the  BASIC  position 
fixing  program  for  the  Commodore  2001.  A  number  of  the 
cursor  control  and  formatting  symbols  on  the  listing  require 
explanation.  These  symbols  are  given  below  by  line  number  at 
first  occurrence. 


130 

Reverse 

field 

'heart'  means  'clear  the  screen' 

130 

Reverse 

field 

1 Q  * 

means 

'move  cursor  down' 

250 

Reverse 

field 

'R' 

means 

'turn  on  field  reverse' 

250 

Reverse 

field 

1  1 

means 

'turn  off  field  reverse' 

250 

Use  of 

the  field 

reverse 

symbols  highlights  the 

'YES'  and  the  'N'  in  'NO'  to  prompt  the  user  to  type  a 
'Y'  or  'N'  on  the  keyboard.  Similar  highlighting  is 
used  on  lines  390,  790,  800  and  810 
740  CHR$ ( 34)  generates  the  double  quote  symbol  to  designate 
'seconds  of  arc' 

740  The  symbol  between  the  double  quotes  which  are  between 
N$(l)  and  N$(2)  is  used  to  designate  ,0',  the  'degrees 
of  arc'  symbol 

910  The  first  symbol  following  the  double  quote  means  'move 
cursor  up',  and  the  next  symbols  each  mean  'move  cursor 
right ' 

The  DATA  statements  on  lines  50001  through  50044  contain 
the  pertinent  parameters  for  each  Loran-C  pair.  The  fields 
contain: 
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1. 


The  Loran-C  station  pair  designator. 

2.  At,  the  sum  of  the  coding  delay  plus  one  way  base¬ 
line  time,  including  the  secondary  phase  correction 
for  an  all  seawater  path,  in  microseconds. 

3.  The  master  station  latitude. 

4.  The  master  station  longitude. 

5.  The  slave  station  latitude. 

6.  The  slave  station  longitude. 

Negative  longitudes  are  West  longitudes.  The  latitudes  and 
longitudes  appear  to  be  in  decimal  form,  but  the  actual  format 
is  DDD.MMSSFF  where 

DDD  designates  degrees, 

MM  designates  minutes, 

SS  designates  seconds, 
and  FF  designates  hundredths  of  seconds. 


50001 

50002 

50003 

50004 

50005 

50006 

50007 
50003 

50009 

50010 

50011 

50012 

50013 

50014 

50015 

50016 

50017 

50018 
50013 

50020 

50021 

50022 

50023 

50024 

50025 

50026 

50027 

50028 

50029 

50030 

5003 1 

50032 

50033 
58034 
50035 
50836 

50037 

50038 
50G39 

50040 

50041 

50042 

50043 

50044 


DfiTfi 
DfiTfi 
IifiTfi 
DfiTfi 
DfiTfi 
DfiTfi 
DfiTfi 
DfiTfi 
DfiTfi 
DfiTfi 
DfiTfi 
DfiTfi 
DfiTfi 
DfiTfi 
DfiTfi 
DfiTfi 
DfiTfi 
DfiTfi 
DfiTfi 
DfiTfi 
DfiTfi 
DfiTfi 
DfiTfi 
DfiTfi 
BfiT  fi 
DfiTfi 
DfiTfi 
DfiTfi 
DfiTfi 
BfiT  fi 
DfiTfi 
DfiTfi 
DfiTfi 
DfiTfi 
DfiTfi 
D  HT  fi 
DfiTfi 
DfiTfi 
DfiTfi 
DfiTfi 
DfiTfi 
DfiTfi 
DfiTfi 
DfiTfi 


4990X, 15972. 23, 16.444395, 
4990V, 34253. 18,16. 444395, 


■169.303120,20. 144916, 
•169. 303 1 20 , 28 . 234 1 77 , 


5930X 

13131. 

So  , 

46 

482720, 

-067 

553 

771. 

41. 

151193 

5930V 

28755. 

02, 

46 

482720, 

-067 

.553 

771. 

46. 

463218 

5990X 

13343. 

, 

51 

575878, 

-122 

220 

224. 

55. 

262085 

5990V 

2892  7  ■ 

36, 

51 

575878, 

-122 

220224. 

47. 

034799 

5990Z 

42266. 

63, 

51 

575878 , 

-122 

220 

224. 

50. 

3&29  r 2 

7930W 

15068. 

02 , 

59 

591727, 

-045 

102 

747. 

64. 

542658 

7930X 

27803. 

f  r  .• 

59 

591727, 

-045 

102 

747. 

62. 

175968 

7930Z 

48212. 

20, 

59 

591727, 

-045 

102 

*!*  A  -J 

<  *+  r  . 

46. 

463218 

7960X 

13804. 

45, 

63 

194281, 

-142 

483 

190. 

,57. 

262021 

7960V 

29651. 

14, 

194281 , 

-142 

4*oo 

190. 

,55. 

262085 

7970W 

30065. 

64 , 

62 

175968, 

-007 

042 

671. 

,54. 

482980 

7970X 

15048. 

10, 

6.i! 

175968, 

-007 

.  042 

671, 

,  68. 

3806 1 5 

7970V 

48944. 

53, 

62 

1 75968 ■ 

-007 

.042 

671 

,64. 

542658 

7970Z 

63216. 

30 , 

62 

175968, 

-007 

,  042 

671, 

.70. 

545261 

7980W , 1 2809 . 54 , 30 . 533874 


7980X 
7980V 
7980Z 
7390 X 
7930V 
7990Z 
S970W 
8370X 
8970V 
993014 
3930X 
9930V 
9930Z 
9940 W 
3940X 
3340V 
9 96 U 14 
3360X 
3960V 
9360Z 
997011 
9370X 
9970V 
99702 
9930X 
3330V 
9390Z 


27443. 33, 
45201 . 88, 
61542.72,30 
12755.97, 

32273. 30, 
50999 . 63 , 

14355.11, 

31 162.06, 
47753.74, 
13635.51, 
36389  •  66.. 

52541 . 3 1 , 
68560. 72, 

13736. 30, 
28094 . 50 , 

41367. 30, 
13797. 20, 

26963 . 93 , 
42221 . 65 • 

57 1 62 . 06 , 

15283.94, 

36685. 12, 

59463. 13, 
30746. 79, 

14875.32, 
32069. 09, 
46590. 10, 


39. 
39 . 


30. 593874 
30. 533874 
33374 
38.522061 
38. 522061, 
38.522061 
10754 
1 0754 
39.510754 
34 . U34604 . 
34 . 034604 . 
34. 034604. 
34 . 084604 

39. 330662 

33. 330662 

39. 330662 
42. 425060 
42. 425060 
42.425060 
42. 425060 
24.48041 
24.48041 
24.43041 
24. 48041 
57 . 090988 
57 . 090388 
57. 030938 


-085 . 1 00330 , 30 . 433302 , 
-085 . 1 00930 , 26 . 3 1 550 1 , 
-035 . 1 00930 ,27.01 5849 , 
-085. 100930, 34.034604, 
016. 430596 , 35.31 2088 , 
016. 430596 , 40 . 582095 , 
016.  4310596 , 42 . 033649 , 
—08 i  .  83 1 2 1 4 ,  80 .  5988  r*  4 , 
-087 .291214 , 42 . 425060 , 
-087 .231214,48. 364984 , 
-077. 544676, 27. 015849, 
-077. 544676, 46. 463218, 
-077 . 544676 ,41.151193, 
-077. 544676, 39.510754, 
- 1 1 8 . 495637 , 47 . 034799 , 
-118. 495637 , 38 . 465699 , 
-113. 495637 , 35 .191818, 
—076. 493386, 46. 482720, 
-076 . 493386 ,41.1511 93 , 
-076.493386, 34. 034604, 
-076. 493386, 39. 510754, 
141.13290  ,24.17077  , 
141.1 9290  ,42. 443700 , 
141.1 9230  , 26 . 362499 , 
141.1 9290  , 09 . 324566 , 
- 170.1 4598 1,52. 494505 , 
-170. 145981,65. 144012, 
-170. 145931 , 57.262021, 


-155. 
-178. 
-069. 
-053. 
-131. 
-119. 
-127. 
-023. 
-007. 
-053. 
-152. 
-131. 
+008. 
+014. 
-023. 
-008. 
-090. 
-097. 
-080. 
-077. 
012. 
027. 
003. 
-085. 
-076. 
-094. 
-080. 
-C53. 
-069. 
-087. 
-119. 
-122. 
-114. 
-067. 
-069. 
—077 . 
-087. 
153. 
143. 
128. 
138. 
+  173. 
-166. 
-152. 


530970 
173020 
583909 
102816 
151965 
443953 
212935 
552175 
042671 
102816 
221122 
151965 
173633 
274700 
552175 
435869 
494360 
500009 
065352 
544676 
312996 
520152 
121590 
100930 
493386 
331847 
065352 
102816 
583909 
291214 
443953 
294453 
481743 
553771 
583909 
544676 
291214 
58515 
430906 
085621 
095523 
105231 
531447 
22 
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